Hyperfine structure of Li and Be 
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A large-scale relativistic configuration-interaction (CI) calculation is performed for the magnetic- 
dipole and the electric-quadrupole hyperfine structure splitting in ^'^Li and ^Be"*". Numerical re- 
sults for the 2^5, 3^5, 2^Pi/2, and 2^ Pi/2 states are reported. The CI calculation based on the 
Dirac-Coulomb-Breit Hamiltonian is supplemented with separate treatments of the QED, nuclear- 
magnetization distribution, recoil, and negative-continuum effects. 
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I. INTRODUCTION 

The hyperfine structure (hfs) of few-electron atoms 
has been an attractive subject of theoretical studies for 
decades, one of the reasons being a few ppm accuracy 
achieved in experiments on Li and Be"*" [lljQ. Despite the 
considerable attention received, a high-precision theoret- 
ical determination of hfs in few-electron atoms remaines 
a difficult task. The main problem lies in the high sin- 
gularity of the hfs interaction and, as a consequence, in 
the dependence of the calculated results on the quality 
of the many-electron wave function near the nucleus. 

The hfs splitting of lithium has traditionally been one 
of the standard test cases for different theoretical meth- 
ods 0. Among various calculations reported in the lit- 
erature, the nonrelativistic ones are the most numer- 
ous; their technique is well developed by now. The best 
numerical accuracy for the nonrelativistic hfs value is 
achieved in variational calculations that use multiple ba- 
sis sets in Hylleraas coordinates 0, d, @ ■ Probably the 
most popular nonrelativistic approach is the multicon- 
figurational Hartree-Fock (MCHF) method 0,[i|, which 
is less computationally intensive but also produces less 
accurate results. The main drawback of the nonrelativis- 
tic methods is that the relativistic effects should be ac- 
counted for separately. There is a way to perform a di- 
rect evaluation of the leading relativistic correction 0, 
but such a calculation is difficult and has not yet been 
done. So far, the relativistic correction was estimated by 
comparing with less accurate calculations based on the 
Dirac-Coulomb-Breit Hamiltonian, or by re-scaling the 
hydrogenic correction. 

There were calculations performed for Li and Be^ with 
the relativistic analog of the MCHF procedure, the multi- 
configurational Dirac-Fock (MCDF) method 0, [U, [11 ■ 
The computational accuracy of these relativistic calcula- 
tions turns out be lower than that of the best nonrela- 
tivistic studies. This is to a large extent due to the fact 
that the electron correlation is more difficult to be ac- 
curately accounted for relativistically than nonrclativis- 
tically. 

Methods that allow a straightforward generalization to 
the relativistic case are many-body perturbation theory 
and its all-order extensions known as the coupled- 



cluster (CC) approach [H, [IE [M [13, [TH . Calculations 
of the Li hfs performed with these methods so far did 
not account for the most part of the triple excitations 
(and, in most cases, for a part of the double excitations 
as well), which led to an incomplete treatment of the 
electron correlation and to a relatively low accuracy of 
the corresponding results. Significant progress in the CC 
calculations was reported in Ref. [ll for the case of Na. 
In that work, all valence triple excitations were included. 
Such an approach, when applied to the Li hfs, would 
significantly improve the accuracy of the CC results. The 
corresponding calculation is presently underway [20| . 

In the current investigation, a relativistic calculation of 
the hfs splitting in Li and Be^ will be performed by em- 
ploying the configuration-interaction (CI) method. Un- 
like the MCDF procedure, the CI method does not in- 
volve a variational minimization and thus is not handi- 
capped by the danger of the variational collapse into the 
negative continuum (which manifests itself in a "sinking" 
of the ground-state energy due to the admixture of the 
negative energy states into the ground-state wave func- 
tion). The CI method has a potential to be more accurate 
than the MCDF method. The only problem is that it 
requires the Dirac spectrum to be sufficiently well repre- 
sented by the model space of one-electron wave functions, 
whereas the MCDF method can produce reasonable re- 
sults with only a few configurations. While this might 
come as a limitation in the case of complicated many- 
electron atoms, the systems at hand, the Li-like atoms, 
are sufficiently simple to be very accurately described by 
the CI method. 

The goal of the present investigation is to perform a 
calculation of the hfs in Li and Be"*" complete to the rel- 



ative order 



where a is the fine-structure constant. 



Such a calculation requires, besides a high-precision de- 
termination of the dominant nonrelativistic contribution, 
a rigorous treatment of the leading relativistic correction 
and the inclusion of the QED effects and ^a"^ . 
Nuclear effects (the recoil and the magnetization distri- 
bution) also contribute on this level. Because of their 
smallness, these effects can be treated nonrelativistically. 
Since we are concerned with the effects of order up to 
only, the Dirac-Coulomb-Breit Hamiltonian may be 
used as a convenient and sound starting point for our 
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investigation. 

The calculation complete to the relative order o? was 
reported for the hfs splitting of the 2^5' and 3^5 states of 
Li and Be+ in our previous paper [^l*!. In this work, we 
extend our calculations to the 2^Pj states and present a 
detailed analysis of various corrections. In particular, the 
recoil correction to the magnetic dipole hfs interaction is 
derived for the case of an arbitrary spin of the nucleus. 
This correction is shown to yield the dominant recoil con- 
tribution for the hfs splitting of the P and higher-Z states 
in medium-Z H-like atoms. To the best of our knowl- 
edge, it has not previously been accounted for in systems 
other than hydrogen and deuterium. 

The paper is organized as follows. In Sec. HH a 
brief summary of the CI method is given. In Sec. IIIH 
we present some basic formulas for the magnetic-dipole 
and the electric-quadrupole interaction and describe our 
CI calculation of the hfs splitting based on the Dirac- 
Coulomb-Breit Hamiltonian. Various corrections to the 
hfs are calculated in Sec. IIVI The results obtained are 
discussed and compared with the experimental data in 
Sec. El 

Relativistic units h — c^ \ and a — (? j (47r) are used 
throughout this paper. 



II. CONFIGURATION INTERACTION 
METHOD 

Relativistic Hamiltonian of an iV-electron atom can be 
written as 



where indices i, j = 1, . . . , iV numerate the electrons, /i^ 
is the one-particle Dirac Hamiltonian, 

hu{i) = cii - Pi + (/3 - 1) m-l- Kuc(n) , (2) 

OL and (5 are the Dirac matrices, V,iuc is the binding po- 
tential of the nucleus, Vci^if) = ck/^'y is the Coulomb 
part of the electron-electron interaction, = jr^ — rj|, 
Vb is the Breit interaction. 



2r,;. 



\a.i ■ OL.j + {OLi ■ -Tij) {OLj ■ fy)] , (3) 



and r = r/r. It is assumed that Hdcb acts in the space 
of the positive-energy eigenfunctions of some one-particle 
Hamiltonian. 

The iV-electron wave function of the system with par- 
ity P, angular momentum quantum number J, and its 
projection M is represented as a linear combination of 
configuration-state functions (CSFs), 



-^{PJM) = Cr^(-frPJM) , 



(4) 



antisymmetrized products of one-electron orbitals ipn of 
the form 



Gn{r) X« 



(r) 



iFn{r) x-K„m,.(r) 



(5) 



where Xftm the spin-angular spinor 22], k = 
(— + 1/2) is the relativistic angular param- 
eter, and m is the angular momentum projection. In the 
present work, we chose the one-electron orbitals V'n to be 
the (positive-energy) eigenfunctions of the one-electron 
Dirac Hamiltonian with the frozen-core Dirac-Fock po- 
tential, 

hDF = Oi-p+{P-l)m + l^„uc(r) + V^p\r) . (6) 

The (non-local) potential V^f^ is defined by its action 
on a wave function, 



V, 



a 
ri2 



[^c{r2)ipiri) -ipc{ri)ip{r2) , (7) 



where the index c runs over the core orbitals. The eigen- 
functions of the Hamiltonian hop form a complete and 
orthogonal basis of one-electron orbitals. 

In the CI method, the ionization energy of the system 
and the mixing coefficients Cr in Eq. ^ are obtained by 
solving the secular equation 



det{{-,rPJM\HDCB\7sPJM) - ErSrs} = 0. 



(8) 



The matrix elements of the Hamiltonian between the 
CSFs can be represented as linear combinations of the 
one- and two-particle radial integrals, 

{jrPJM\ Hdcb\1sPJM) = ^ drs{ah) I{ah) 

ah 

+ a ^ ^ {ahcd) [i?^ {ahcd) + i?f [ahcd)] . 



k ahcd 



(9) 



Here, a, fe, c, and d specify the one-electron orbitals, d^s 
and v'^rf are the angular coefficients, I{ab) are the one- 
electron radial integrals, and R'^{ahcd) and R^{abcd) are 
the Coulomb and Breit two-electron radial integrals. The 
radial integrals are defined by 



nab) , 



(10) 



{ab\Vc,B\<^d) = a ^ 



knik 
k rrik 

'ia nio Jc -nic ^ ]d md,]b -rat " 



X C^.Z': . _™ C^?„\ Rf^iabcd) , 

(11) 



where 7^ denotes the set of additional quantum numbers 
that determine the CSF. The CSFs are constructed as 



where C:? ™ , ^ are the Clebsch-Gordan coefficients. 
After integrating over angular variables, the expression 
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for the one-particle integral reads 
d K 



I{ab) 



dr 



dr r J \ dr r 



2mFaFb 



(12) 



The Coulomb integral is given by 



oo k 
' < 



X J dridr2-j^Wacin)Wbd{r2), (13) 

where Wab — GaGb + FaFb and C^'^^ is the spherical ten- 
sor with components c[j\r) = \J Am j (2J -|- 1) Yjuij")- 
The expression for the Breit integral is more complex; 
it can be found in Ref. [1^. The angular coefficients 

dra and v^) can be evaluated analytically 0, [1^. In 
the general case, formulas for them are rather cumber- 
some. A number of packages is available in the literature 
for the numerical evaluation of the angular coefficients 



III. HYPERFINE SPLITTING 

A. Magnetic dipole hyperfine splitting 

The relativistic Fermi-Breit operator of the magnetic 
dipole hyperfine interaction is given by 



47r 



(14) 



where /x is the operator of the nuclear magnetic moment, 
acting in the nuclear subspace. The operator T^^^^ acts 
in the electronic subspace; it is given by the sum of the 
one-electron operators t'^^^{i), 



(15) 



In the nonrelativistic limit, the operator t^^^ turns into 



''NR' 



1 

m 



3f(s ■ f) — s 8n ^, . 

— + Y^W^ 



(16) 



where I and s are the one-electron operators of orbital 
angular momentum and spin, respectively. The three 
terms in the brackets are often referred to as the orbital, 
the spin-dipole, and the Fermi contact term, respectively: 



''NR — ''I 



.(1) 



(17) 



The relativistic value of the energy shift due to the 
magnetic dipole hyperfine interaction is obtained as the 
expectation value of the Fermi-Breit operator on the wave 
function of the system with atomic angular momentum 



F and its projection Mp. Employing the Wigner-Eckart 
theorem, the nuclear variables are separated and inte- 
grated out and the energy shift is represented in terms 
of the reduced matrix element of the operator T^^\ 



AEmi = {FMf\Hmi\FMf) = 



|e| jj_ 
An 21 

X [F(F+l)-/(/+l) - J(J+1)] 



v/J(J+l)(2J+l) 



(18) 



where /i is the magnetic moment of the nucleus, fi — 
{II\ij,q\II) , I is the nuclear spin, and J is the total an- 
gular momentum of the electron. Experimental data for 
the magnetic dipole hfs splitting are usually expressed 
in terms of the hyperfine interaction constant Aj, which 
does not depend on F, 



A 



AE 



Ml 



\e\ M (J||tW||J) 



{FMf\I ■ J\FMf) 47r / ^ J( J + 1)(2 J + 1) ' 

(19) 

The reduced matrix element of the operator T^"^' 
should be evaluated with the CI many-electron wave 
functions (|3]), obtained by solving the secular equation 
([8]). Matrix elements of the operator T^^^ between indi- 
vidual CSFs can be expressed as linear combinations of 
matrix elements of the one-electron operator t^^^ between 
the single-particle orbitals, 

(7.PJ||T(i)||7.FJ) = '^rK<^b) (a||t(^)||6) , (20) 

where a and b numerate the one-electron orbitals and d'^rs 
are the angular recoupling coefficients. Packages for the 
numerical evaluation of the coefficients d^'r) are available 
in the literature [2^, |2^ . The reduced matrix element of 
the one-electron operator t^^^ is given by 

(a||*(l)||&) ^ -{na + Kb){-Ka\\C^^^\\Kb) 



drr-^{GaFb + FaGb). (21) 



The energy shift due to the hfs splitting can easily be 
calculated for the hydrogen-like ion. In this case, the sum 
([20| consists of a single term and the radial integral in 
Eq. (PT|) is calculated analytically (in the point-nucleus 
limit). In the present work, we will need the nonrela- 
tivistic limit of Eq. ()2ip for the hydrogen-like ion, which 
reads 



2{Zaym^ 1 2j + 1 



n3 2l + iyj{j + l) 



(22) 

Using this result, it is convenient to introduce the follow- 
ing parametrization of the magnetic hyperfine constant 
Aj: 



A-, 



1 



TOp ^jv /J(J + 1)(2L+ 1) 



Gmi{Z) , 



(23) 
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where n is the principal quantum number of the valence 
electron, nip is the proton mass, and /ijv — |e|/(2TOp) 
is the nuclear magneton. The function Gmi{Z) is di- 
mensionless; its numerical value is unity for a hydrogen- 
like nonrelativistic atom in the point-nucleus and non- 
recoil limit. Gaii is a slowly varying function of the nu- 
clear charge number Z and the quantum numbers J and 
L, which is convenient for the representation of numer- 
ical results. This definition of the function Gmi differs 
slightly from the one used in our previous work |2l| by 
the fact that it does not include the nonrelativistic mass 
scaling factor {l + m/M)~^ . We presently choose to treat 
this part of the recoil effect (also referred to as the nor- 
mal mass shift) on an equal footing with the other recoil 
corrections. A parametrization similar to that in Eq. (j23p 
was previously used in Refs. [30i, i31i]. 

Numerical results of nonrelativistic calculations are of- 
ten presented in terms of the orbital (a;), the spin-dipole 
(asd), and the Fermi contact (oc) hyperfine parameters, 
induced by the three terms in Eq. (fTT]) and defined as Q 



ai ^ {LSMlMs\Y,^-^\LSMlMs) , (24) 

i=l 

asd = {LSMM J2 ° ^ ^'> \LSMlMs) , 

(25) 

a, = {LSMlMs\ J2 ^^S!i^i^\LSMLMs) , (26) 



i=l 



(27) 



with Ml = L and Ms = S. The connection of the (non- 
relativistic limit of the) function Gmi with the hyperfine 
parameters expressed in atomic units is given by 



Gmi{Z) 



n^J{J +1){2L+1) 



2Z3 

The coefficients Ci are 



(c/ ai + Csd asd + Cc ttc) ■ 
(28) 



ci 



Csd = 



(L-J) 



LJ{J +1) ' 

3(5 ■£) {L-J) - L{L + 1){S ■ J) 
5L(2L- 1)J(J+ 1) 
(S-J) 



3SJ{J+1) 



(29) 
(30) 
(31) 



where 



(L-J) 
(S-J) 
(S-L) 



[J{J + 1) + L{L + 1) - S{S + l)]/2 , (32) 
[J( J + 1) - L{L + 1) + S{S + l)]/2 , (33) 
[J( J + 1) - L{L + 1) - S{S + l)]/2 . (34) 



B. Electric quadrupole hyperfine splitting 

The scalar part of the interaction between an electron 
and the nucleus is given by 



^{'^7 '^Pl J ■ ■ ■ I ) 



1 



(35) 



where r and are the coordinates of the electron and 
the jth proton, respectively. Averaging this interaction 
over the internal nuclear coordinates and using the stan- 
dard multipole expansion of |r — Vp^ \ ~^, one obtains (see 
Ref. [H for the details) 



CO „ 

Kv(r-,e,$) EE (F)i„tcrn = / 



dr' r'' 



Piir') 



,1+1 



r'-r) + —9{r - r') 



(36) 



where (• • •) denotes the averaging, 8 and (f> are the angles 
that fix the orientation of the intrinsic nuclear system 
with respect to the laboratory frame, and the nuclear 
charge density component pi is defined as 



dfp{r)c'i\r) 



(37) 



The first term in Eq. (|36|) {I — 0) yields the stan- 
dard Coulomb interaction between the electron and the 
nucleus with an extended charge distribution. The term 
with I — 1 vanishes after averaging with the electron wave 
function of a definite parity. The term with I = 2 gives 
rise to a splitting of the energy level (of an electronic 
state with J > 1/2), known as the electric quadrupole 
one. The corresponding interaction is conveniently writ- 
ten in the form 



H£2 = aT(2).g(2). 



(38) 



(2) 

Here, Qav is the operator of the nuclear quadrupole mo- 
ment averaged over the internal (radial) nuclear coordi- 
nates, 



/ intern 



ArC(2)(e,$). (39) 



The normalization constant N is 
rR 



N= / drr''p2ir), 



(40) 







where R is the nuclear radius. The operator T*^^-' acts on 
the electronic variables. It is given by 

T(2)=^<(2)(z) = -5:/(r.)C(2)(n), (41) 
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where the radial distribution function /(r) is 



1 

, r > R, 



fir) 



1 

N 



dr'r' p2{r')^, r<R. 



(42) 



The distribution function /(r) can easily be calculated 
analytically for several simple models of the nuclear- 
charge distribution. So, if p2 does not depend on r within 
the nucleus, P2(?') 9{R — r), 



/(0 = ;^(l + 51n|) , r<R. (43) 



If P2{r) (X S{R — r), then 



/(r) = ^, r<R. 



(44) 



In the point-quadrupolc limit, the function f{r) takes the 
standard form, /(r) = r^^. 

The finite nuclear size effect is very small for the elec- 
tric quadrupole splitting and its inclusion in calculations 
is not necessary at present. However, we observed that 
the usage of the extended charge distribution consider- 
ably improves the stability and the convergence of nu- 
merical calculations. The reason for this is that the ex- 
tended distribution removes the singularity of the 
point-quadrupole interaction. 

Using the standard technique of the angular- 
momentum algebra (see, e.g., Ref. Q), the energy shift 
due to the electric quadrupole interaction can be ex- 
pressed in terms of the reduced matrix elements of the 
electronic operator T^^^. The correction to the energy 
is usually expressed in terms of the hyperfine structure 
constant Bj, which docs not depend on the total angular 
momentum of the system F, 



_ |C(C + l)-/(/-H)J(J + l) ^ 
^^^^ 2/(2/-l)J(J+l) ' 

where C = F{F + 1) - /(/ + 1) - J(J + 1) and 



(45) 



Bj = 2Q , 



2J(2J- 1) 



(2J+ l)(2J + 2)(2J + 3) 



(j||r(^)||j) 



(46) 



Here, Q is the nuclear quadrupole moment, defined as 



g = (/A/|^(3z, 



rl)\IM)M=i = 2{II\Q^^^\II) 



(47) 

The reduced matrix element of the operator T^^^ 
should be evaluated with the many-electron wave func- 
tions ^ , obtained by the solution of the secular equation 
([S]). Matrix elements of the operator T^^^ between indi- 
vidual CSFs can be expressed as linear combinations of 



matrix elements of the one-electron operator i^^^ between 
the single-particle orbitals, 

J| |T(2) I \j,PJ) ^ J2 di'J («l 1*^'^ 1 1^> , (48) 

a<b 



where a and b numerate the one-electron orbitals and dl-s 
are the angular recoupling coefhcients [2G] . The reduced 
matrix element of the one-electron operator t^^^ is given 
by 

(a||t(2)||6) = -(kJ|C(2)||«;,) / drf{r){GaGb+FaFh). 

Jo 

(49) 

Similarly to the magnetic dipole hyperfine constant 
Aj, the electric quadrupole hyperfine constant Bj can 
be conveniently parameterized by introducing the dimen- 
sionless function Ge2, which turns into unity for a nonrel- 
ativistic hydrogen-like ion in the point-nucleus and non- 
recoil limit. 



B.J = 



a{Zaym-^ 2 J - 1 



1 



i3 J+1 k(k + 1)(2L + 1) 



GE2iZ). 



(50) 

where k = {-iy'+^+^/^{J + 1/2). 

Results of nonrclativistic calculations are often ex- 
pressed in terms of the quadrupole parameter 6^, defined 
as 

bg = (LSMlMsI J2 ° 3 "' \LSMlMs)m,=l.Ms=s ■ 

i=l 

(51) 

The connection between the (nonrelativistic limit of the) 
function Ge2 and the parameter bq expressed in atomic 
units is given by 



Ge2{Z) 



TT'iJ + 1)k{k+1){2L + 1) 



Z3(2J- 1) 

where the coefficient Cq is Q 

_ 6(X • J)2 - 3(i • J) - 2L{L + 1) J(J + 1) 
^'^ ~ L(2i- 1)(2J + 3)(J+ 1) 

C. Details of the CI calculation 



i-Cqbq), (52) 



(53) 



To perform a CI calculation, we devised a code, in- 
corporating and adapting a number of existing packages 
[i4, 26, 27|, m, 1^ for setting up the CSFs, calculating 
angular-momentum coefficients, and diagonalizing the 
Hamiltonian matrix. The largest number of CSFs simul- 
taneously handled was about a half of a million. A careful 
optimization of the code was necessary to keep the time 
and memory consumption of the calculation within rea- 
sonable limits. Care was taken to prevent re-calculating 
the angular-momentum coefficients for the pairs of CSFs 
that differ by the principal quantum number of a single 
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electron only. An optimized ordering of CSFs allows one 
to drastically reduce the number of angular-momentum 
coefficients to be evaluated. A similar optimization was 
introduced in the calculation of the Coulomb and Breit 
radial integrals. The radial integrals with the same pair 
of electron states in the innermost radial integration were 
grouped together and evaluated simultaneously. 

The dominant part of the hfs splitting is delivered by 
the Dirac-Coulomb Hamiltonian. This is the most de- 
manding part of the calculation, since a high relative 
precision is required. One of the factors defining the ac- 
curacy of the calculation is the quality and the size of the 
space of one-electron orbitals from which the CSFs are 
constructed. We take this space to be a part of the finite 
basis set of eigenvectors of the Dirac equation, obtained 
by the dual-kinetic-balance method [34] and constructed 
with _B-splines [ssj . 

For a given number of B-splines Ua, all eigenstates 
were taken with the energy < e < mc^{l + ZaE-^^^ 
and the orbital quantum number I < Z^ax, where 
the value of -Emax was varied between 0.5 and 6 and 
^max, between 1 and 7. Three main sets of one- 
electron orbitals were employed in the present work: 
(A) 20s 20p 19d 19/ 18g 18h with Ua = 44 and S^ax = 
3.0, (B) 14sl4pl4dl3/133l3/il2zl2fc with = 34 
and iSniax = 0.5, and (C) 25s25p24d with = 54 
and Smax — 6.0. Here, the notation, e.g., 20p means 
20pi/2 20p3/2. Calculational results were first obtained 
with the set (A) and then corrected for contributions of 
the higher partial waves with the set (B) and for a more 
complete representation of the Dirac spectrum with the 
set (C). The computation became rather intensive for the 
P states, so the basis set (A) was reduced to include the 
states with I < 3 only in this case. Usage of several sets of 
one-electron orbitals allowed us to efficiently control the 
completeness of the representation of the Dirac spectrum 
in our calculations. 

The analysis of the convergence of the partial- wave ex- 
pansion was performed by identifying increments of the 
results induced by the increasing cutoff parameter Zmax- 
The omitted tail of the expansion was estimated by a 
polynomial least-square fitting of the increments in l/l. 
In most cases, the error due to the termination of the 
expansion was found to yield the largest uncertainty to 
the Dirac-Coulomb hfs value. 

The set of the CSFs employed in the calculation was 
obtained by taking all single, double, and triple excita- 
tions of the reference configuration with at least one elec- 
tron orbital with I < 1 present. The contribution of the 
remaining triple excitations was found to be negligible 
for the S states. For the P states, it was estimated by 
repeating the calculation with a smaller basis but with 
the above restriction replaced hy I < 2. 

Inclusion of the Breit interaction into the Dirac- 



Coulomb Hamiltonian yields only a small correction in 
the case of Li and Be+. Because of this, it is sufficient 
to use much smaller basis sets for its evaluation, which 
simplifies the computation greatly. The Breit-interaction 
correction was obtained as the difference of the CI results 
with and without the Breit interaction included into the 
Hamiltonian, evaluated with the same set of CSFs. 

Results of our CI calculations of the magnetic dipole 
and the electric quadrupole hfs splitting are presented in 
TablesHlandim respectively. The CI values obtained with 
the Dirac-Coulomb Hamiltonian are listed under the en- 
try "Coulomb" ; the entry "Breit" contains the correction 
due to the inclusion of the Breit interaction. The com- 
parison presented in the tables demonstrates significant 
deviations of our CI values from the MCDF results by 
Bieroh et al. [lol. [ll[ and from the CC results of John- 
son et al. jl8j|. In the case of Be"*", Ref. [llj reports 
estimations of the calculational errors, considered by the 
authors to be the conservative ones, but our CI results 
are well out of these error bars for all the states studied. 

The deviation from the MCDF calculations is the 
strongest for the quadrupole splitting in Be"^. In this 
case, our CI value differs from the MCDF one already 
in the second digit, while the claimed accuracy of the 
MCDF result is about 10~^. A similar deviation is ob- 
served also for the quadrupole splitting in Li. At the 
same time, agreement of our calculations with the non- 
relativistic studies [1, [s^ is much better, on the level 
of lO^'^. This observation leads us to a conclusion that 
the MCDF results for the quadrupole splitting are, most 
probably, in error. A possible explanation for this is that 
the highly singular point-quadrupole interaction might 
lead to considerable numerical errors when evaluated on 
approximate relativistic wave functions. In our calcu- 
lations, we detected such problems; they were solved by 
using the extended charge distribution for the quadrupole 
interaction. 

In order to make possible a detailed comparison with 
high-precision nonrelativistic results available in the lit- 
erature, we have to identify the nonrelativistic part of 
our CI values. This was achieved by repeating the full 
set of the CI calculations for different values of the fine- 
structure constant a (namely, three values with ratios 
a' /a = 0.9, 1, and 1.1 were used). For each value of a, 
the finite nuclear-charge distribution correction was eval- 
uated (as described in the next section) and subtracted 
from the CI values. The point- nucleus results thus ob- 
tained were fitted to a polynomial in a, assuming the 
absence of the linear term. In this way, the CI results 
with the physical value of a were separated into three 
parts: the nonrelativistic point-nucleus contribution, the 
relativistic correction, and the finite nuclear-charge cor- 
rection. 
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TABLE I: Individual contributions to the magnetic dipole hfs splitting in ^'^Li and ^Be"*", in terms of the function Gjv/i if not 
specified otherwise. For ®Li, only the contributions different from those for ^Li are listed. The values of the nuclear magnetic 
moments are taken from Ref. |33]- The entries are labelled as follows: "NR(point)" denotes the point-nucleus nonrelativistic 
result; "Relativistic" is the total relativistic correction; "Coulomb" is the relativistic hfs value obtained with the Dirac-Coulomb 
Hamiltonian; "Breit" is the Breit-interaction correction; "BW" is the nuclear magnetization distribution correction; "NMS" is 
the normal mass shift; "SMS" is the specific mass shift; "SO" is the sum of the normal and specific spin-orbital recoil corrections 
induced by Eqs. (|6H) and (|62|l : "Negative-energy" is the contribution of the negative-energy part of the Dirac spectrum. 







2^S 


3^5' 






Ref. 






J_j1 










NR(point) 




0.215251 


0.168340 


0.073905 


-0.024348 








0.215254 (4) 


0.168351 (13) 






Hylleraas [5] 






n 91 ^1 Q 




u.u * ooy 






J. l_,dcl \j L V ±Oul\j 




000205 


000159 


000018 


—0 000036 


Coulomb 




0.215385 (5) 


0.168440 (9) 


0.073923 (1) 


—0.024364 (10) 








0.21527 




0.07396 


-0.02476 


MCDF [10] 






0.21565 


0.16861 


0.07389 


-0.02425 


CCSD \18\ 


Breit 




0.000016 


0.000016 


-0.000003 


0.000000 




QED 




0.000182 (4) 


0.000143 (3) 


0.000048 (1) 


-0.000085 (2) 




BW 




—0.000024 (5) 


-0.000019 (4) 


—0.000002 


0.000009 (2) 




Recoil 


NMS 


-0.000050 


-0.000039 


-0.000017 


0.000006 






SMS 


0.000002 


0.000002 


0.000027 


—0.000055 






SO 


0.000000 


0.000000 


-0.000001 (1) 


-0.000002 (1) 




Negative-energy 




0.000002(1) 


0.000002 (1) 


-0.000003(1) 


-0.000003 (2) 




Total 




0.215512 (8) 


0.168544 (11) 


0.073972 (2) 


-0.024493 (10) 




Total" fMHzl 

L J 




401.755(15) 


93.095 (6) 


45.966(1) 


-3.044(1) 




Experiment [MHz] 




401.7520433(5)'' 


93.106(11)'' 


45.914(25)'* 


-3.055 (14)'' 












46.010(25)"= 














46.024(3)-'' 










6 Li 










Coulomb 




0.215382 (5) 


0.168438 (9) 


0.073922 (1) 


-0.024363 (10) 




BW 




-0.000022 (13) 


—0.000017 (10) 


-0.000002 (1) 


0.000008 (5) 




Recoil 


NMS 


—0 000059 


— 000046 


—0 000020 


000007 






SMS 


0.000002 


0.000002 


0.000031 


-0.000064 




Total 




0.215504(14) 


0.168538(14) 


0.073974 (2) 


-0.024502(11) 




TotaF [MHz] 




152.122(10) 


35.250(3) 


17.4058 (5) 


-1.1530 (5) 




Experiment [MHz] 




152.136839(2)'' 


35.263(15)= 


17.375 (18)'' 


-1.155 (8)'' 












17.386 (31^ 














17.394(4)-'' 




















NR(point) 




0.390544 


0.335066 


0.221132 


0.00989 








390549 (Q\ 








±± y iiCl CLCLB 1 1 






39050 


33504 


0.22113 


00967 


MCHF |8] 


Relativistic 




0.000664 


0.000563 


0.000162 


-0.00015 


Coulomb 




0.391030 (6) 


0.335468 (9) 


0.221302 (2) 


0.009800 (25) 








0.39094 (4) 




0.22140(1) 


0.0091 (4) 


MCDF [11] 


Breit 




0.000039 


0.000042 


-0.000021 


-0.000001 




QED 




0.000289 (12) 


0.000248 (10) 


0.000137(5) 


-0.000181 (7) 




BW 




-0.000062 (6) 


-0.000053(5) 


-0.000005 (1) 


0.000027 (3) 




Recoil 


NMS 


-0.000071 


-0.000061 


-0.000040 


-0.000002 






SMS 


0.000002 


0.000002 


0.000057 


-0.000080 






SO 


0.000000 


0.000000 


-0.000007(2) 


-0.000018 (5) 




Negat ive- energy 




0.000005 (3) 


0.000005 (2) 


-0.000009 (4) 


-0.000011 (6) 




Total 




0.391233(15) 


0.335651 (15) 


0.221413(7) 


0.009533 (27) 




TotaP [MHz] 




-625.08 (2) 


-158.897(7) 


-117.919 (4) 


-1.015(3) 




TotaP [MHz] 




-625.11 (3) 


-158.905(7) 


-117.925 (4) 


-1.016(3) 




Experiment [MHz] 




-625.00883705 (1)* 




-118.6(36)' 







" nCLi) = 3.2564268 (17) 
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Beckmann et al, 1974 [Ij, 
" Bushaw et al, 2003 [38], 

Orth et al, 1975 [39], 

Walls et al, 2003 [40], 
^ Das and Natarajan, 2008 [41], 
f /xC'Li) = 0.8220473 (6), 

Orth et al, 1974 [42], 
^ MC'Be) = -1.177432 (3), 
^ M(^Be) = -1.177492 (17), 

Wineland et al, 1983 [2], 
' Bollinger al, 1985 [43]. 



For the P states and the magnetic dipole hfs, the non- 
relativistic limit of the CI results needs to be separated 
into three parts, corresponding to the three terms of the 
nonrelativistic decomposition of the hfs operator (jl7p . 
To this end, we carried out identical calculations both 
for the relativistic magnetic dipole hfs operator ([15]) and 
for the spin-dipole and the orbital parts of its nonrel- 
ativistic decomposition. Applying the fitting procedure 
described above, we identify the nonrelativistic limit of 
the CI values as well as the spin-dipole and orbital hfs 
parameters a; and a ad- The remaining contact parame- 
ter Oc is then unambiguously deduced. (We prefer not to 
perform a direct calculation for the contact term since the 
corresponding operator contains a (5-function and needs a 
regularization when evaluated on relativistic wave func- 
tions.) 

The nonrelativistic hfs parameters obtained in this way 
are listed in Table [1111 The nonrelativistic results for the 
2^P state were obtained from the relativistic calculations 
for the 2^Pi/2 state. Since in the present work the hfs pa- 
rameters are needed for the purpose of comparison only, 
we do not assign the uncertainty to them (which is diffi- 
cult to do reliably since they are obtained by a fit). The 
cornparison with the previous nonrelativistic calculations 
@) H) HE] presented in the table exhibits a remarkably 
good agreement of our values with the high-precision re- 
sults obtained in a Hylleraas-type calculation by Yan et 
al [|. 

In Tables HI and [m the entry "NR(point)" labels the 
nonrelativistic, point-nucleus limits of the functions Gmi 
and Ge2 obtained by the fitting procedure described 
above. Because of the fitting, the uncertainties are not 
ascribed to them; we expect that they are somewhat less 
accurate than the corresponding relativistic values. The 
comparison is drawn with the most accurate previous 
nonrelativistic calculations. A much better agreement 
is observed with the previous nonrelativistic results than 
with the relativistic ones. 



IV. CORRECTIONS TO THE HYPERFINE 
SPLITTING 

While the evaluation of the relativistic hfs value is the 
most computationally intensive part of the calculation, a 
high-precision theoretical determination of the hfs split- 



ting requires inclusion of a number of important correc- 
tions. In this section, we present a detailed description 
of each of them in turn. 



A. QED effects 

For the magnetic dipole hfs splitting, the leading (in 
Za) QED contribution originates from the anomalous 
magnetic moment of the electron g^. The effect is ac- 
counted for by multiplying the spin-dependent terms in 
Eq. mi) by ge/2 w a/{2T:), see, e.g., Ref. [sj. So, the 
leading QED correction to the function Gmi is given by 



27r 



GMl,sd{Z) + Gmi,c{Z) 



(54) 



where GMi,sd and Gmi.c are the contributions to the 
function Gmi induced by the spin-dipole and the contact 
term in Eq. (|16p . respectively. 

The higher-order terms of the Za expansion (the bind- 
ing QED corrections) induce important contributions 
and should be taken into account alongside with the lead- 
ing effect. The binding corrections to the contact term 
can be written in a form analogous to that for the hydro- 
gen hfs Hi, El, 



6G 



QED.bind 
Ml 



(Z) = ^ Gmi AZ) I ZaiT (ln2 - ^ 



+ {Zaf 



8 .2 



In [Za) + ln(ZQ;) -|- 



(55) 



The coefficients 021 and 020 are different from the hydro- 
genic case and not known at present. One can, however, 
use their hydrogenic values as crude estimates. In our 
calculations, we will use the results for the hydrogenic 2s 
state, a2i = —1.1675, 020 = 11.3522 0,14^, and assume 
a 100% uncertainty for them. This treatment of the QED 
effects coincides with those of Refs. [1, 0] but is different 
from other previous investigations, where the binding ef- 
fects were continually neglected. Such neglect can hardly 
be justified since the higher-order terms change the total 
QED contribution by 40% for lithium and by 60% for 
beryllium. 

The binding corrections to the spin-dipole and orbital 
parts of hfs are relevant for the states with I > only. 
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TABLE II: Individual contributions to the electric quadrupole hfs splitting of the "i^P-^j^ state, in terms of the function Ge2 if 
not specified otherwise. The notations are the same as in Table U The values of the nuclear quadrupole moments are taken 
from Ref. ;37]. 

^Li Ha ^Be+ Rcf\ 



INK 


U.0o02d0 


0.050260 


O.i 1 IlW 






0.0498 




0.1717 


MCHF [8] 




0.0498 




0.1727 


FCPC [36] 


Relativistic 


-0.000004 


-0.000004 


-0.000013 


Coulomb 


0.050260 (3) 


0.050260 (3) 


0.172150 (7) 






0.051085 




0.18356(3) 


MCDF [10, nj 


Breit 


-0.000004 


-0.000004 


-0.000024 


QED 


0.000000 (2) 


0.000000 (2) 


0.000000(11) 




NMS 


-0.000012 


-0.000014 


-0.000031 




SMS 


0.000012 


0.000014 


0.000030 




Total 


0.050256 (4) 


0.050256 (4) 
-0.0044(1)'' 


0.172125 (13) 




Total [MHz] 


-0.216 (4)" 


2.281 (16)'= 




Experiment [MHz] 


-0.221 (29) 


-0.010(14) 




[39,42] 


" QCIa) = -40.55 (80) mbarn, 



Q(®Li) = -0.82 (2) mbarn, 
" (5(^Be) = 52.88 (38) mbarn, 



TABLE HI: Non-relativistic hfs parameters, in a.u. 



Ion 


2^5 


3^S 




22p 






Ref. 












0,1 






Li 


2.90589 


0.67336 


-0.21467 


-0.013477 


0.063125 


-0.022617 


This work 




2.90592(5) 


0.67341(5) 


-0.21478(5) 








Hylleraas [5] 




2.9051 


0.6731 


-0.2151 


-0.01346 


0.06311 


-0.02239 


MCHF [8] 




2.903 


0.6745 


-0.2136 


-0.01341 


0.06309 


-0.02242 


FCPC [36] 


Be+ 


12.4974 


3.1769 


-1.0842 


-0.10269 


0.48520 


-0.18362 


This work 




12.4976(3) 












Hylleraas [5] 




12.496 


3.1767 


-1.0856 


-0.10265 


0.48516 


-0.18310 


MCHF [8] 




12.493 


3.181 


-1.070 


-0.1020 


0.4851 


-0.1842 


FCPC [36] 



They enter in the relative order q;(Zq;)^ and are presently 
unknown. Numerical calculations for the hydrogenic case 
[4^ show that their nominal order can be enhanced by 
the second power of logarithm. We thus estimate the 
uncertainty due to their neglect by 

\GMi.sd + Gmi,i \ — {ZaY \n{Za) . 



For the electric quadrupole splitting, the QED correc- 
tion has not been calculated so far. Its relative nominal 
order is a{Za)^. According to our analysis, this correc- 
tion diverges in the point-quadrupole limit i? 0, which 
means that the nominal order is enhanced by Ini? « 6. 
We, therefore, estimate the error in Ge2 due to the ne- 
glect of the QED effects by multiplying it by the factor 
of 

lQa{Zaf . 



B. Nuclear recoil 

Within the nonrelativistic approach, the nuclear re- 
coil effect on the energy levels and on the wave functions 
of the system is accounted for by introducing two addi- 
tional terms in the Hamiltonian, traditionally referred to 
as the normal mass shift (NMS) and the specific mass 
shift (SMS). They are given by 

i 

respectively, where M is the mass of the nucleus. Alter- 
ations of the wave function due to the additions to the 
Hamiltonian give rise to the corresponding corrections to 
the hfs. The NMS part of the recoil can be factorized 
out and expressed in terms of the reduced mass. It is ac- 
counted for by multiplying the nonrelativistic hfs value 
by a factor of (1 + m/M)~^ [13] . The inclusion of Hnms 
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into the CI Hamiltonian leads to the same effect but adds 
some relativistic corrections. The SMS part of the recoil 
effect is to be evaluated numerically, by incorporating 
HsMS into the CI Hamiltonian and by identifying the 
corresponding alteration of the hfs splitting. 

It should be stressed that, despite the fact that our 
original CI Hamiltonian is the relativistic one, the inclu- 
sion of the operators Hnms and Hsms in it does not 
fully account for the relativistic recoil effects, since the 
operators themselves are obtained within the nonrela- 
tivistic approximation only. This fact was often disre- 
garded in the past, e.g., in Ref. |4S]. The complete treat- 
ment of the leading [~(Za)^] relativistic recoil correction 
to energy levels of the system is achieved by employing 
the operator [49l | 

ij 

(58) 

Numerical calculations with this operator were per- 
formed, e.g., in Ref. In our present investigation, 
the relativistic recoil effects are negligible as compared 
to other sources of the theoretical uncertainty. We thus 
use the nonrelativistic operators for the description of the 
recoil effects. 

Matrix elements of Hsms between the individual CSFs 
can be expressed in terms of the angular coefficients 
introduced in Eq. with the multipolarity k = 1, 

{jrPJM\HsMshsPJM) = 

X 4^'{abcd) V{ac) V{bd) , 

abed 

(59) 

where the radial integrals are (see, e.g., Ref. [isj ) 

V{ac) = (Ka||CW||Kc) / dr 
Jo 

Ka{Ka + 1) - K^K,, + 1)" 



d 
dr 



2r 



Fa 



d_ _ Ka{Ka ~ 1) - Kc{Kc - 1) 

dr 2r 



Fc } . (60) 



So far, we discussed the recoil corrections to hfs that 
are induced by the wave functions. There are, however, 
also recoil corrections to the hyperfine interaction itself. 
The recoil correction to the magnetic dipole hfs interac- 
tion arises through the spin-orbit coupling in the scalar 
component of the nuclear current. This correction de- 
pends on the spin of the nucleus /. In the case of hydro- 
gen {I — 1/2), it was derived many years ago in Ref. [Slj, 
whereas in Ref. [sl] it was reported for the case of deu- 
terium (/ = 1). To the best of our knowledge, this cor- 
rection was previously unknown for the arbitrary spin of 
the nucleus and was not accounted for in calculations of 
the hfs of systems other than hydrogen and deuterium. 



The spin-orbital (SO) recoil correction to the magnetic 
dipole hfs interaction is obtained in Appendix by using 
the expression for the current of a particle with an ar- 
bitrary spin derived in Ref. [53| . The result is repre- 
sented by Eq. (|A.7[) . It can be conveniently split into the 
normal (SON) and specific (SOS) parts, analogously to 
the normal and the specific mass shift of energy levels, 
F[so — HsoN + HsoS: with 



ff^^o^-2^(5-l)J-E^^' (61) 

i * 

Hsos = ^ (5 - 1) I • E *sosi^, j) , (62) 

W(*,j) = ^-3^ + ^-3^- (63) 



Here, / is the operator of the nuclear spin, Zj is the op- 
erator of the orbital angular momentum of ith electron, 
and g is the g- factor of the nucleus, 



^l M 1 

/^AT TTlp I 



(64) 



The SON interaction is proportional to the orbital 
part of the nonrelativistic decomposition of the magnetic 
dipole hfs operator. It is easy to see that this part of the 
SO recoil effect can be accounted for by modifying the 
orbital hyperfine parameter a/ by 



ai 



ai 



TO „ q — 1 

1 + — Z- 

M g 



(65) 



It is interesting to note that, comparing to the NMS 
effect, the SON correction is enhanced by a factor of Z, 
which makes it a dominant recoil effect in the hfs splitting 
of medium-Z H-like ions (for electronic states with Z > 0). 
For Li2+, the ratio of the SON and the NMS effects is 
—0.5 for the 2pi/2 state and —1.2 for the 2^3/2 state. For 
the lithium-like systems, however, the SON and the SOS 
corrections tend to cancel each other, the net effect being 
rather small numerically. 

In the present work, we calculate the SOS correction 
by using perturbation theory to the lowest order. For 
the electronic configuration with a single valence elec- 
tron beyond the closed core shell, the contribution to the 
function Gmi due to the SOS effect can be expressed as 



'Ml 



v?J{J +1){2L 
Zm g — 1 



1) 



E(-I)(c^^l^sosol«c), (66) 



where v denotes the valence electron state with the an- 
gular momentum projection /i^, = 1/2, c is the core elec- 
tron state with the angular momentum projection /x^, 
and tsoSo is the zeroth spherical component of the oper- 
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ator tsos- The radial integral is evaluated to yield 



where 



is the Zemach moment obtained by fold- 



{cv\tsoso\'"c) 



6 



/ + l)(2j, + 1) 
111 

Jv Jv Jc 



U{cv) V{vc) , 



(67) 



where 



U{cv) = {k,\\C^'^\k,) / drr-^G,G, + F,F,), 
Jo 

(68) 

and V{vc) is defined by Eq. (|60p . It is easy to see that 
the radial integral (p7|) vanishes for the S states. 

The calculational results for the individual recoil con- 
tributions to the magnetic dipole hfs splitting are listed 
in Table U under entries "NMS", "SMS", and "SO". The 
results obtained for the NMS and SMS parts are in good 
agreement with the previous evaluations of these correc- 
tions. The entry "SO" represents the sum of the SON 
and the SOS corrections. Because of a large cancella- 
tion between these two parts, we calculate both of them 
by perturbation theory. The uncertainty specified in the 
table was evaluated by comparing results obtained with 
different potentials in the zeroth-order Hamiltonian. 

The scalar component of the nuclear current yields 
also a correction to the electric quadrupole interaction 
[s^ . [5^. This correction is induced by the nuclear 
spin and can be interpreted as a shift of the nuclear 
quadrupole moment (see Appendix for details). The in- 
duced contribution is included into the observable value 
of the nuclear quadrupole moment and thus is not needed 
to be taken into account in the theoretical description of 
the electric quadrupole hfs. 



C. Nuclear size and magnetization distribution 

Due to a high singularity of the hfs interaction at the 
origin, the nuclear structure effects (particularly, the dis- 
tribution of the nuclear magnetic moment) have signifi- 
cant influence on the magnetic dipole hfs and should be 
taken into account in atomic calculations. An accurate 
theoretical description of these effects is a demanding 
problem. A way for its rigorous solution was paved in 
recent studies [13, [H, HI] ■ Practical realizations of this 
approach, however, are so far restricted to two- and three- 
nucleon systems j54 . i55i] and their extension to more com- 
plex nuclei like ^Li and ^Be looks problematic. 

The most widely used approach up to now is to ac- 
count for the extended nuclear magnetization distribu- 
tion [the Bohr-Weisskopf (BW) effect] by means of the 
Zemach formula (STj . According to the original formula- 
tion, the nuclear correction to the magnetic dipole hfs of 
an S state of an H-like atom is represented by a simple 
multiplicative factor. 



ing together the electric charge Pe{r) and magnetization 
Pmir) densities 



= J drdr' pe{r) Pm{r')\r - r'\ 



(70) 



Formula accounts for both the charge and the 
magnetization distribution. Since the charge distribution 
effect is usually taken into account in a more complete 
way by modifying the Coulomb nuclear potential in the 
Hamiltonian it should be subtracted from the total 
Zemach correction. The finite nuclear charge (FNC) cor- 
rection is obtained from Eq. (|69|) by setting {r)em — (^)e, 
where (r)e is the electric charge radius defined as 



drpe{r) 



(71) 



More detailed studies of the FNC correction in H-like 
atoms with including the relativistic effects were reported 
in Refs. ^]. For lithium and beryllium, the relativis- 
tic effects are small and enter mainly through the alter- 
ation of the exponent in the Za and (r)e dependence by 
terms ~ (Za)'^. 

Using the hydrogenic result [s^] for the exponent of the 
Za and (r)e dependence, we write the generalization of 
the nonrelativistic FNC correction in the form valid for 
an arbitrary state of few-electron atoms 



SG 



IK^{Z) = -2 {Za {r),r'-^ GmiAZ) , 



(72) 



{Z)^-2Za (r)e™ Gmi(^), 



(69) 



where 7 = \Jl — [ZaY and Gmi.c is the contact part of 
Gmi- Using Eq. ((72)) . one should keep in mind that the 
charge radius (r) e is different from the charge root-mean- 
square (rms) radius (r^)y^, which is usually listed in ta- 
bles. The conversion factor depends somewhat on the 
model of the nuclear charge distribution. For the Gaus- 
sian model, the connection is (r)e = -^/S/ (Stt) (r^)^^. 
The values of the rms radii of the nuclei ^'^Li and ^Be 
were taken from Ref. |59] . The rms radii and the numer- 
ical results for the FNS correction obtained by Eq. ((72|) 
are presented in Table IIVI (the relative values of the cor- 
rection, SG^jl*^ /Gmi.c, are hsted). An independent eval- 
uation of the FNS correction was performed, by repeating 
the CI calculations for different values of the nuclear ra- 
dius and by fitting the increments to the analytical form 
([7^ . The FNS correction obtained in this way agreed 
very well with the analytical results presented in the ta- 
ble. 

In the present investigation, separate values of the FNS 
correction are not necessary since this effect is already 
included in the CI part of the calculation. However, we 
use the values of the FNS correction in order to identify 
the point-nucleus limit of our results (particularly, for 
the comparison with the point-nucleus results of Ref. [1]) 
and for improving the stability of the fit in extracting the 
nonrelativistic limit of our calculations (the FNC correc- 
tion is the only part of the CI values that is linear in 
a). 
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There is no need to specify explicitly to which elec- 
tronic states the results for the nuclear corrections in 
Table IIVI correspond, because the relative values of the 
corrections are listed. Our numerical calculations show 
that, with a good accuracy, the relative values of the nu- 
clear corrections do not depend on the particular state. 
(Of course, for the P states, the relative value should be 
evaluated with respect to the contact part of the correc- 
tion.) 

The Zemach correction induced by the magnetization 
distribution (the BW effect) can be written in a form 
valid for an arbitrary state as 



-2Za [(r),„. - (r),] Gmi,c(^) ■ (73) 



The Zemach radius is usually not tabulated and should 
be derived from data available for the charge and mag- 
netization rms radii according to Eq. (|70p . with an addi- 
tional input of the distribution models. For the Gaussian 
model 



P{r) = Pa exp(-Ar^) 



(74) 



employed for the charge and magnetization distributions, 
the Zemach radius is readily obtained analytically, 



,1/2 



(75) 



For more sophisticated distribution models, one has to 
evaluate radial integrations in Eq. (j70p numerically. In 
order to test the model dependence of the Zemach radius 
(with fixed values of the charge and magnetization rms 
radii), we performed its numerical evaluation with the 
two-parameter Fermi model. The same results as for the 
Gaussian model are obtained, which leads us to conclude 
that the model dependence is negligible. 

The values listed in Table IIVI for the magnetic rms ra- 
dius are the average of data tabulated in Ref. [B^I and 
the errors are their mean-square deviation. Under the 
entry "Zemach", we tabulate the numerical results for 
the BW correction obtained by Eq. ([75)1 : the error as- 
cribed to them originates from the uncertainties of the 
magnetization and charge radii. 

The second approach to the description of the BW ef- 
fect considered in the present work is based on the single 
particle (SP) model of the nuclear magnetic moment and 
will be referred to as the SP approach in the following. 
Within the SP model, the nuclear magnetic moment is 
assumed to be induced by the odd nucleon (proton, when 
Z and A are odd and neutron, when Z is even and A is 
odd) . The odd nucleon is assumed to have an effective g 
factor, which is fixed so that it yields the experimental 
value of the nuclear magnetic moment. The treatment 
of the magnetization distribution effect on hfs within the 
SP model was originally developed in Refs. [6l|, [g^l and 
later in Ref. [13] ■ The spin-orbit interaction of the odd 



nucleon was introduced into this approach in Ref. [Slj. 
Our present treatment closely follows the procedure de- 
scribed in Refs. [sil. [63l|. 

The wave function of the odd nucleon is assumed to 
satisfy the Schrodinger equation with the central poten- 
tial of the Woods-Saxon form and the spin-orbital term 
included (see, e.g., Ref. (U) 

Vir) = -VoHr) + —(bsoir) l-cr + Vcir), (76) 



where 



4>soir) 



Vso dT{r) 



Anipr dr 



nr) = 



exp 



r-R 



(77) 



(78) 



and Vc is the Coulomb part of the interaction (absent 
for neutron), with the uniform distribution of the charge 
{Z — 1) over the nuclear sphere. The parameters Vq, Vso, 
R, and a were taken from Ref. (6^ . where they were ob- 
tained by fitting electron scattering data. The fitting was 
not evaluated for ^Be, so we use the parameters for its 
closest odd-neutron neighbour, ^^C. For ^Li, the nuclear 
spin is integer (/ = 1), and so one needs to make an addi- 
tional assumption about the value of the orbital angular 
momentum of the nucleus. We used the value L = [6^1 
in our calculations. 

The nuclear magnetic moment can be evaluated within 
the SP model to yield [IH 



1 



:9s 



I 



1 



2/+ 1 



9L 



PN 



2(/+l) 



9s 



2 4(/+l 

7(2/ + 3) 2/+1 



2(/+l) 4(/+l) 



for / = L + - 



9L ■ 



for I = L - 

2 ' 
(79) 

where / and L are the total and the orbital angular mo- 
mentum of the nucleus, respectively, is the g factor 
associated with the orbital motion of the nucleon {gL ~ 
1 for proton and gL = for neutron) and gs is the effec- 
tive nucleon g factor, determined by the condition that 
Eq. (1751) yields the experimental value of the magnetic 
moment. 

It was demonstrated in Ref. [3l| that, within the SP 
model, the BW effect can be accounted for by adding a 
multiplicative magnetization-distribution function F{r) 
to the standard point-dipole hfs interaction psp . The 
distribution function is given by j63j 
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TABLE IV: Nuclear parameters (in fm) and the relative values of the nuclear corrections (in ppm) to the magnetic dipole hfs. 
The abbreviations are as follows: "FNS" denotes the nuclear charge distribution correction, "Zemach" labels the BW correction 
obtained with the Zemach formula, "SP" labels the BW correction evaluated within the SP approach, "const" indicates that 
the odd-nucleon wave function is taken to be a constant within the nucleus. 



Isotope 






{r)em — {r)e 


FNS 




Bohr-Weisskopf effect 




Total nuclear 




[fm] 


[fm] 


[fm] 




Zemach 


SP SP (const) 


Total 


correction 


«Li 


2.54(3) 


3.12(22) 


1.38(19) 


-268(3) 


-160(20) 


-50 -89 


-100(60) 


-368(60) 


'Li 


2.43(3) 


2.80(8) 


1.19(8) 


-257(3) 


-135(9) 


-112 -99 


-112(23) 


-369(23) 


^Be 


2.52(1) 


2.67(6) 


1.06(5) 


-356(4) 


-160(8) 


-158 -161 


-158(16) 


-514(16) 



F{r) = ^ Tdr' 
A* Jo 



■gs 



I - 



1 21+1 



2 4(/+l) 



MAT 



21-1 
'8(/+l) 



9s 



1 21+1 

2 ^ 4(/+l) 



for / = i + 1/2 and 

F(r) = ^ r dr'r'^\u{r')\^ 
M Jo 



2(1+1) 



9s 



dr' r'^|u(r')|^ — 



2/ + 3 
8(7+1) 



9s 



1(21 + 3) 21+1 2^ , M 

WTYy-WTT)'^'''^'^^^''^ 

1(21 + 3) 21+1 



2(1+1) 4(/+l) 



r (t)so(r) 9L 



(80) 



(81) 



for / = i — 1/2. In the above formulas, u(r) is the wave 
function of the odd nucleon. It can easily be seen that 
F(r) = 1 outside the nucleus. 

In the present work, we evaluated the BW correction 
within the SP approach as described above. In addition, 
we considered a simplified version of this approach ob- 
tained by assuming the wave function of the odd nucleon 
to be just a constant within the nucleus. By compar- 
ing the two corresponding results, we can conservatively 
estimate the dependence of the SP values on the param- 
eters employed in the Woods-Saxon potential. Calcula- 
tional results for the BW correction obtained within the 
SP approach are listed in Table HVl under the labels "SP" 
(the full SP approach) and "SP(const)" (the SP approach 
with the constant wave function of the odd nucleon). 

The total results for the BW correction listed in the 
table were obtained as follows. For ^Li and ^Be, we em- 
ploy the results of the SP model as the final values. The 
uncertainty for ''Li was taken to be the largest deviation 
from the final result. For ^Be, all three values fall very 
close to each other, so we assign the 10% uncertainty to 
the final result. The nucleus ^Li has an odd neutron and 
an odd proton; one thus can hardly expect it to be de- 
scribed well by the SP model. In this case, we use the 
plain average of the three values; the error was chosen so 
that it covers all three results. 

We would like to stress that there are nontrivial nu- 
clear structure effects, which are ignored both within the 
SP model and within the Zemach approach. Since these 
effects cannot be reasonably estimated at present, our 



uncertainties of the BW correction yield the order of the 
expected error only. 

The final values of the BW correction and their un- 
certainties are listed in Table |T] under the entry "BW" . 
Our results for ''Li and ^Be are reasonably close to the 
Zemach- formula values of Ref . Q . The result of Ref. [llj 
for ^Be is larger than ours by a factor of 4. This is 
because the authors of Ref. [IJj] used (r)m instead of 
{r)em - {r)e m Eq. In Ref. 0, the BW correction 

was evaluated within the SP approach with the constant 
odd-nucleon wave function. The corresponding results 
for ^'^Li nearly coincide with our values obtained within 
the same approach. In Refs. [H, [6^, the BW correction 
was calculated by using the Fermi [l^j or the uniform 
[66| distribution of the magnetization density over the 
nucleus. Their results fall between our values obtained 
with different models. In most of other previous studies, 
the BW effect was not accounted for. 



D. Negative continuum 

The negative-continuum (NC) contribution might be 
of some importance in calculations involving the opera- 
tors that mix the upper and the lower components of the 
Dirac wave function. The magnetic dipole hfs operator 
is of this kind, so we have to obtain an estimation for 
this correction. In the present investigation, we calculate 
the NC contribution by employing perturbation theory 
to the first order. 
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For the electronic configuration with a single valence 
electron beyond the closed core shell and to first order in 
the electron-electron interaction, the NC correction can 
be written as 



n3J(J+ 1)(2L+ 1) 



-1 £„<0 

2E 



E 



[{vc\VcB\nc) - {cv\VcB\nc)] (n|4^^|^ 

£"7. — Eji 



E 



[{vc\VcB\vn) - {cv\VcB\vn)] (n|4'^|c) 



{v\U\n) {n\t^^'\v) 



(82) 



where v denotes the valence electron state with the angu- 
lar momentum projection — 1/2, c is the core electron 

state with the angular momentum projection /ic, io^"* is 
the spherical component of the magnetic dipole hfs oper- 
ator defined by Eq. US]), Vcb ^Vc + Vb is the sum 
of the Coulomb and Breit parts of the electron-electron 
interaction, and the summation over n is performed over 
the negative-energy part of the Dirac spectrum. The 
states V, c, and n are assumed to be eigenvectors of a 
single-particle Dirac Hamiltonian h with the screening 
potential C/, h — a. ■ p + {(3 — l)m + Vnuc(?') + U{r). 

We mention that the NC correction may depend 
strongly on the choice of the Hamiltonian h. It is, there- 
fore, important to use the same single-particle Hamilto- 
nian for the evaluation of the NC correction as in the 
CI part of the calculation. So, in the present investiga- 
tion the screening potential U in Eq. ([5^ was fixed as 

'-' ^ ^DF ■ 

It should be also noted that, evaluating Eq. (|5^ . one 
cannot neglect the Breit part of the electron-electron in- 
teraction as compared to the Coulomb one. For the 
negative-energy part of the Dirac spectrum, both of these 
interactions induce contributions of the same order of 
magnitude. 

Formula ([82)1 for the NC correction ignores contribu- 
tions of the second and higher orders in the electron- 
electron interaction. Their unambiguous description is 
possible within QED only. For lithium and beryllium, 
the electron correlation is strong and the perturbation ex- 
pansion converges slowly. We thus assign the uncertainty 
of 50% to the NC contribution obtained by Eq. 



V. DISCUSSION 

The calculational results for the magnetic-dipole hfs 
splitting of the S^S*, 2'^Pi/2, and 2^^3/2 states of 

^'^Li and ^Be^ are listed in Table HI expressed in terms 
of the dimensionless function Gmi defined by Eq. (^5)1 . 
The values presented for the S states differ from those 
in our previous work ^2ll| in two ways. First, we now 



treat the normal mass shift as a correction, rather than 
by including it into the definition of the function Gmi- 
(Of course, this difference has no effect on the total the- 
oretical prediction for the hfs or the hyperfine constant 
Aj.) Second, wc perform a more detailed analysis of the 
BW effect, and so the uncertainty of this correction is 
changed. 

In order to convert the function Gmi into the hyper- 
fine constant Aj, an additional experimental input in the 
form of the magnetic moment of the nucleus is needed. 
This issue might contain some ambiguities since the tab- 
ulated values of the nuclear magnetic moments [13] are 
often inconsistent. In the case of ^Li and ^Li, we as- 
sume the values originally obtained by Beckmann et al. 

to be the most reliable ones. For ^Be, the choice is 
less obvious, and we present the theoretical results for 
the hyperfine constant Aj for two different experimental 
values of the magnetic moment. 

For '^Li, we observe good agreement of our theoretical 
predictions with all the experimental results listed, ex- 
cept the recent measurement of the 2^Pi/2 state by Das 
and Natarajan |4l|, which claims to be accurate to better 
than 0.01%. The theoretical prediction is away from this 
measurement by about IScr and we presently see no way 
to explain this deviation theoretically. 

Our theoretical prediction for the ground-state hfs 
splitting of ^Li is in slight disagreement with the high- 
precision experimental result. For the ground state of 
beryllium, the deviation is larger and amounts to 3 or 
4 a, depending on the value of the nuclear magnetic mo- 
ment used. There are two possible explanations of these 
discrepancies: underestimated systematic effects in the 
experimental values of the nuclear magnetic moments 
and nontrivial nuclear-structure effects in the theoretical 
predictions. Basing on the experimental data available, 
we cannot unambiguously distinguish between these two 
explanations. It would have been possible if the hfs split- 
ting of two different states were accurately measured for 
the same isotope. The existing measurements of the hfs 
of excited states, however, are not yet sensitive to the in- 
consistencies in values of the nuclear magnetic moment. 

The comparison presented for the hfs of excited states 
of lithium and beryllium indicates that our theoretical 
predictions are more accurate than the experimental re- 
sults. The general agreement with the experimental data 
for the excited states is good, the only exception being 
the results of Ref. [4l|. 

Our calculational results for the electric quadrupole hfs 
splitting of the 2^P3/2 state in ^'^Li and ^Be+ are listed in 
Table [m expressed in terms of the dimensionless function 
Ge2 defined by Eq. (|50|) . It is remarkable that all theo- 
retical contributions to the electric quadrupole splitting 
seem to be well under control, so that the resulting theo- 
retical predictions for the function Ge2 are obtained with 
very good accuracy. If accurate experimental investiga- 
tions of the quadrupole splitting were possible for Li-like 
ions, they would lead to a high-precision determination 
of the nuclear quadrupole moments, which are difhcult to 



15 



measure directly. In the absence of such investigations, 
the theoretical predictions for the hyperfine constant Bj 
are obtained by using the tabulated values of the nuclear 
quadrupole moments [37|. Our results are in agreement 
with the scarce experimental data available. 
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APPENDIX: CORRECTIONS TO THE 
HYPERFINE SPLITTING DUE TO THE SCALAR 
COMPONENT OF THE NUCLEAR CURRENT 

Let us first consider the simplest case of the electron 
and the nucleus being the spin-1/2 Dirac particles. The 
electron-nucleus interaction is then given by the standard 
scattering amplitude [gJ 



(A.l) 



where D^'' is the photon propagator and je and j„uc 
are the electromagnetic current of the electron and the 
nucleus, respectively, 

= ^u{p') [(P + Pr + ^^''"lA u{p) . (A.2) 

Here, u is the free Dirac spinor, m is the mass of the 
particle, q = p' — p, and cr^'' = («/2)[7^, 7^]. Expressing 
the time component of the nuclear current in terms of 
the free spinors w in the rest frame, one arrives at [6^ 



Jnuc = w* ( 1 - 



8A/2 



icr ■ q X p 



(A.3) 



where M is the nuclear mass. The first term in the brack- 
ets in the above expression corresponds to the standard 
Coulomb interaction between electron and the nucleus. 
The third term represents the spin-orbital coupling and 
induces a recoil correction to the hfs we are interested 
here. 

The generalization of the expression (jA.2p for the case 
of the nucleus with an arbitrary spin was obtained in 
Ref. [5^. The time component of the current reads 



1 

2M 



i^ip') [Fe (E + E') + G„ 7°r • q] , 

(A.4) 



where and Gm are the electric and magnetic form 
factors of the nucleus, respectively. 



S 

-S 



(A.5) 



and the vector S is constructed from components Si, 
which are generalizations of the Pauli matrices. After 
expressing the current in terms of the spinors ^0 in the 
rest frame, one gets the generalization of Eq. (|A.3p to the 
case of an arbitrary-spin nucleus [53| . 



7° 

Jnuc 



' SO 



Fe — {2Gm — Fe) 



.I-{qxp) 



2AP 



(A.6) 



where / is the operator of the spin of the nucleus. 
The form factors are identified as [53]: Fe{0) — 1 and 
Grn{0) ~ 5/2, where g is the g factor of the nucleus. 

The third term in the brackets of Eq. (|A.6p induces 
a first-order (in the electron-nucleus mass ratio) recoil 
correction to the magnetic dipole hfs splitting. Taking 
into account that, in the center-of-mass system, the total 
momentum of the atom is zero and transforming this 
term into the coordinate space (see Ref. [13]), we obtain 
the interaction of the form 



Hso — 



Za 
2M2 



n X Pi 



(A.7) 

where indices i and j numerate the electrons in the atom. 
For the hydrogen atom, / = (1/2)(t, Eq. (|A.7p repro- 
duces the well-known result of Ref. A similar re- 
coil correction to the Zeeman splitting of multi-electron 
atoms was reported in Ref. [68j . 

The second term in the brackets of Eq. (|A.6p can be 
split into the contact and the quadrupole part, which in- 
duce corrections to the Lamb shift and to the quadrupole 
hfs splitting, respectively. Both of these corrections were 
evaluated in Ref. [5^. The result for the quadrupole in- 
teraction due to the second term in Eq. (IA.6 ) is (with 
the additional factor of 2, corrected in Ref. 



SH 



E2 



Za{g 



2M2 



■i^AV,V,i 



3 ^ 



where 



A 



1/(2/— 1) , I is integer, 
1 / (2/) , / is half-integer . 



(A., 



(A.9) 



Using the standard angular-momentum algebra, we 
transform Eq. (jA.Sp into the form analogous to Eq. 



5HE2--aT^'^^ ■5Q^'^\ (A.IO) 

where JQ^^-* is the correction to the operator of the 
quadrupole moment. 



5Q 



(2) 



^(g-l)A 
2M2 



\/6(/(g)/)(2) . 



(A.11) 
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Taking into account Eq. (|47)) . the correction to the nu- 
clear quadrupole moment is identified, which is 



5Q 



Z{g-l)I 

M2 ' 

Z(g-l)(/-l/2) 
A'P 



I is integer , 



/ is half-integer . 



(A.12) 



The numerical values of the induced quadrupole mo- 
ment for the isotopes considered in this work are: 
<5Qind(^Li) = —0.39 mbarn, (5Qind(^Li) = —0.15 mbarn, 
and 5(3ind(^Be) = —0.18 mbarn, to be compared with 
the total values of the nuclear quadrupole moments ■ 
QCU) = -40.55 (80) mbarn, Q{^U) = -0.82 (2) mbarn, 



Q(9Be) = 52.88 (38) mbarn. 

It should be mentioned that the correction (jA.12[) does 
not have immediate experimental consequences. It van- 
ishes for the nuclear spin / = and 1/2. Nuclei with 
the spin / > 1/2 have a quadrupole moment, and so the 
correction (jA.12p appears only together with the "pure" 
nuclear quadrupole moment. If the values of the nuclear 
quadrupole moments are derived from experimental ob- 
servations, the induced correction is included in them 
and thus does not have to be accounted for in theoretical 
descriptions of the quadrupole splitting. It should be in- 
cluded, however, when the nuclear quadrupole moments 
are calculated basing on microscopic nuclear models, as, 
e.g., in Ref. [t^- 
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